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ABSTRACT 

We explore the effects of dust in cosmologically distributed intervening galaxies on 
the high redshift universe using a generalised model where dust content evolves with 
cosmic time. The absorbing galaxies are modelled as exponential disks which form 
coevally, maintain a constant space density and evolve in dust content at a rate that 
is uniform throughout. We find that the inclusion of moderate to moderately weak 
amounts of evolution consistent with other studies can reduce the mean observed B- 
band optical depth to redshifts z>\ by at least 60% relative to non-evolving models. 
Our predictions imply that intervening galactic dust is unlikely to bias the optical 
counts of quasars at high redshifts and their evolution in space density derived there- 
from. 
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1 INTRODUCTION 

The recent discovery of large numbers of quasars at ra- 
dio and X-ray frequencies with very red optical-to-near- 
infrared continua suggests that existing optical surveys may 
be severely incomplete (eg. Webster et al. 1995 and refer- 
ences therein). Webster et al. (1995) and Masci (1998) have 
argued that the anomalous colours are due to extinction by 
dust, although the location of the dust remains a highly con- 
troversial issue. Intervening dusty galaxies which happen to 
lie along the line-of-sight of otherwise normal blue quasars 
are expected to redden the observed optical continuum, or if 
the optical depth is high enough, to remove quasars from an 
optical flux-limited sample (eg. Wright 1990). As suggested 
by existing obervational and theoretical studies of cosmic 
chemical evolution however (Pei & Fall 1995 and references 
therein), one expects a reduction in the amount of dust to 
high redshift. Consequently, one then also expects that the 
probability of a background object being either reddened or 
obscured to be reduced. 

The effects of foreground dust on observations of objects 
at cosmological distances has been discussed by Ostriker & 
Heisler (1984); Heisler & Ostriker (1988); Fall & Pei (1989, 
1993); Wright (1986, 1990) and Masci & Webster (1995). Us- 
ing models of dusty galactic disks, these studies show that 
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the line-of-sight to a high redshift quasar has a high proba- 
bility of being intercepted by a galactic disk, particularly if 
the dust distribution is larger than the optical radius of the 
galaxy. Based on the dust properties of local galaxies, it is 
estimated that up to 80% of bright quasars to z ~ 3 may be 
obscured by dusty intervening systems. The principle issue 
in these calculations was that realistic dust distributions in 
galaxies which are 'soft' around the edges, will cause many 
quasars to appear reddened without removing them from a 
flux-limited sample. 

None of the above studies however considered the ef- 
fects of evolution in dust content. Cosmic evolution in dust 
is indirectly suggested by numerous claims of reduced chem- 
ical enrichment at z>2. Evidence is provided by observa- 
tions of trace metals and their relative abundances in QSO 
absorption- line systems to z ~ 3 (Meyer & Roth 1990; 
Savaglio, D'Odorico & Moller 1994; Pettini et al. 1994; Wolfe 
et al. 1994; Pettini et al. 1997; Songaila 1997), which are 
thought to arise from intervening clouds or the haloes and 
disks of galaxies. These studies indicate mean metallicities 
~ 10% and ^1% solar at z ~ 2 and z ~ 3 respectively, 
and dust-to-gas ratios ^8% of the galactic value at z ~ 2. 
These estimates are consistent with simple global evolution 
models of star formation and gas consumption rates in the 
universe (Pei & Fall 1995). If the observed metallicities in 
QSO absorption systems are common, then their interpreta- 
tion as galactic disks implies that substantial evolution has 
taken place since z ~ 3. If the quantity of dust on cosmic 
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scales also follows such a trend, then one may expect the 
effects of obscuration to high redshift to be reduced relative 
to non-evolving predictions. 

In this paper, we continue to model the effects of inter- 
vening galactic dust on the background universe at optical 
wavelengths using a more generalised model where the dust 
content evolves. We explore the effects of our predictions on 
quasar number counts in the optical and their implication 
for quasar evolution. 

This paper is organised as follows: The next section 
briefly describes the generalised model and assumptions. 
Section ^ describes the model parameters and their values 
assumed in our calculations. Model results are presented and 
analysed in Section ^. Implications on quasar statistics and 
evolution are discussed in Section |^. Other implications are 
discussed in Section [] and all results are summarised in Sec- 
tion 0. Unless otherwise stated, all calculations assume a 
Friedmann cosmology with qo = 0.5, and Hubble parameter 
/150 = 1 where Ho = 50/iso kms" 1 Mpc - . 



2 THE EVOLUTIONARY DUST MODEL 

We calculate the probability distribution in total dust opti- 
cal depth from model galaxies along any random line-of-sight 
as a function of redshift by following the method presented 
in Masci & Webster (1995). This was based on a method in- 
troduced by Wright (1986) which did not include any effects 
of evolution with redshift. Here we generalise this model by 
considering the possibility of evolution in the dust proper- 
ties of galaxies. In the discussion below and unless otherwise 
indicated by a subscript, we define r to be the total optical 
depth encountered by emitted photons and measured in an 
observer's B bandpass (effectively at A = 4400 A). 

We assume the following properties for individual ab- 
sorbing galaxies. Following previous studies (eg. Wright 
1986, Heisler & Ostriker 1988), we model galaxies as ran- 
domly tilted exponential disks, where the optical depth 
through a face-on disk decreases exponentially with distance 
r from the center: 



~(r,z) = 7-0(2) e r/r °. 



(1) 



ro is a characteristic radius and to (2), the value of r through 
the center of the galaxy (r = 0). The redshift dependence 
of to is due to the increase in absorber rest frame frequency 
with redshift. 

Since we wish to model the observed _B-band optical 
depth to 2<6, we require an extinction law £(A) = t\/tb 
that extends to wavelengths of ~ 630A . We use the analyt- 
ical fit for £(A) as derived by Pei (1992) for diffuse galactic 
dust in the range 500A<A<25/im. The optical depth in an 
observer's frame through an absorber at redshift 2 (to (2) in 
equation Fl) can be written: 



T (2) = T S £ 



Ar 



1 + 2 



(2) 



where tb is the rest frame B-band optical depth through 
the center of an individual galactic absorber. 



2.1 Evolution 

Equation (Q) must be modified if the dust content in each 
galaxy is assumed to evolve with cosmic time. The opti- 
cal depth seen through the center of a single absorber at 
some redshift, to (2), will depend on the quantity of dust 
formed from past stellar processes. For simplicity, we assume 
all galaxies form simultaneously, maintain a constant space 
density, and increase in dust content at a rate that is uni- 
form throughout. We also assume no evolution in the dust 
law £(A) with redshift. Even though a lower mean metal- 
licity at high redshift may suggest a different wavelength 
dependence for the dust law, there is no evidence from lo- 
cal observations of the diffuse ISM to support this view (eg. 
Whittet 1992). 

We parameterise evolution in dust content by follow- 
ing simulations of the formation of heavy metals in the cold 
dark matter scenario of galaxy formation by Blain & Longair 
(1993a, 1993b). These authors assume that galaxies form by 
the coalescence of gaseous protoclouds through hierarchical 
clustering as prescribed by Press & Schechter (1974). A fixed 
fraction of the mass involved in each merger event is con- 
verted into stars, leading to the formation of heavy metals 
and dust. It was assumed that the energy liberated through 
stellar radiation was absorbed by dust and re-radiated into 
the far-infrared. They found that such radiation can con- 
tribute substantially to the far-infrared background intensity 
from which they use to constrain a model for the formation 
of heavy metals as a function of cosmic time. Their models 
show that the comoving density of heavy metals created by 
some redshift 2, given that star formation commenced at 
some epoch zsf follows the form 



fi m (z) oc In 



ZSF 



where 2 < zsf- 



(3) 



We assume that a fixed fraction of heavy metals con- 
dense into dust grains so that the comoving density in dust, 
£ld(z), follows a similar dependence as equation (^). The 
density in dust relative to the present closure density in no 
exponential disks per unit comoving volume is given by 



a-; 



n M d 

Pc 



(4) 



where p c = 3Ho/8ttG and Md is the dust mass in a single 
exponential disk. This mass can be estimated using Eq.7-24 
from Spitzer (1978) where the total density in dust, pd, is 
related to the extinction Av along a path length L in kpc 
by 



(Pd) 



1.3 x 10" 



PB ( J±f ) iAv/L). 



(5) 



p g and e are the density and dielectric constant of a typical 
dust grain respectively and the numerical factor has dimen- 
sions of gmcm -2 - see Spitzer (1978). Using the exponential 
profile (equation |l|) where r(r) oc Av(r) and integrating 
along cylinders, the dust mass in a single exponential disk 
can be found in terms of the model parameters tb and ro. 
We find that the comoving density in dust at some redshift 
scales as 



Q. d {z) oc TB(z)nor , 



(6) 



where tb(z) is the central B-band optical depth and ro 
the dust scale radius of each disk. Thus, the central opti- 
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Figure 1. Optical depth in an observer's B-bandpass as a func- 
tion of redshift through a single model absorber defined by equa- 
tion (ph. rg is the rest frame central B-band optical depth and 
Zdust t ne dust formation epoch. 



cal depth, tb{z), in any model absorber at some redshift is 
directly proportional to the mass density in dust or heavy 
metals as specified by equation (^): 



tb(z) oc In 



1 + zsf 
1 + z 



(7) 



The redshift dependence of optical depth observed in 
the fixed B-bandpass due to a single absorber now involves 
two factors: first, the extinction properties of the dust as 
defined by equation (^) and second, its evolution specified 
by equation (hj). The star formation epoch zsf can also be 
interpreted as the redshift at which dust forms. From here 
on, we therefore refer to this parameter as Zdust - a hypoth- 
esised "dust formation epoch" . By convolving equations (^|) 
and (^), and requiring that locally: tq(z = 0) = tb, the 
observed optical depth through a single absorber at some 
redshift z < Zdust now takes the form: 



T (z) 



TB 



<(3 



+ Z 



ln(l + z) 



(8) 



ln(l + Zdust 

Figure ^ illustrates the combined effects of evolution and 
increase in observed frame B-band extinction with redshift 
defined by equation (^|). The extinction initially increases 
with z due to a decrease in corresponding rest frame wave- 
length. Depending on the value for Zdust, it then decreases 
due to evolution in dust content. This latter effect dominates 
towards z dust . 

The characteristic galactic dust radius ro defined in 
equation (|l| is also given a redshift dependence in the sense 
that galaxies had smaller dust-haloes at earlier epochs. The 
following evolutionary form is adopted: 



ro(z) = r (1 + z) 



5 < 0, 



(9) 



where 5 gives the rate of evolution and ro is now a 'local' 
scale radius. Evolution in radial dust extent is suggested by 
dynamical models of star formation in an initially formed 
protogalaxy (Edmunds 1990 and references therein). These 
studies show that the star formation rate and hence metal- 
licity in disk galaxies has a radial dependence that decreases 
outwards at all times. It is thus quite plausible that galaxies 



have an evolving effective 'dust radius' which follows chem- 
ical enrichment from stellar processes. 

Our parameterisation for evolution in galactic dust 
(equations ^ and |^) is qualitatively similar to the 'accre- 
tion models' for chemical evolution of Wang (1991), where 
the effects of grain destruction by supernovae and grain for- 
mation in molecular clouds is taken into account. The above 
model is also consistent with empirical age-metallicity rela- 
tionships inferred from spectral observations in the Galaxy 
(Wheeler, Snedin & Truran 1989), and models of chemical 
evolution on a cosmic scale implied by absorption-line ob- 
servations of quasars (Lanzetta et al. 1995; Pei & Fall 1995). 



3 MODEL PARAMETERS AND 
ASSUMPTIONS 

3.1 Model Parameters 

Our model depends on four independent parameters which 
describe the characteristics and evolutionary properties of 
intervening galaxies. The parameters defined 'locally' are: 
the comoving number density of galaxies no, the character- 
istic dust radius ro, and dust opacity rg at the center of an 
individual absorber. The evolution in tb and ro is defined by 
equations (Q) and ^ respectively. Parameters defining their 
evolution are S for ro, and the 'dust formation epoch' Zdust 
for tb- Both no and ro have been conveniently combined 
into the parameter r 9 where 

T g = n 7rr — , (10) 

with being the Hubble length. This parameter is pro- 
portional to the number of galaxies and mean optical depth 
introduced along the line-of-sight (see Section^). It also rep- 
resents a 'local' covering factor in dusty galactic disks - the 
fraction of sky at the observer covered in absorbers. 

In all calculations, we assume a fixed value for no. From 
equation (|l^) , any evolution in the comoving number density 
no is included in the evolution parameter S for ro (equation 
Thus in general, 8 represents an effective evolution pa- 
rameter for both ro and no. Our model is therefore specified 
by four parameters: r s , tb, 5 and z du st- 



3.2 Assumed Parameter Values 

Our calculations assume a combination of values for the pa- 
rameters (t 9 , tb) and (8, z dua t) that bracket the range con- 
sistent with existing observations. The values (r g , tb) are 
chosen from previous studies of dust distributions and ex- 
tinction in nearby spirals. From the studies of Giovanelli 
et al. (1994) and Disney & Phillipps (1995) (see also refer- 
ences therein) we assume the range in central optical depths: 
0.5^tb<4, while dust scale radii of 5jS(ro/kpc)<30 are as- 
sumed from Zaritsky (1994) and Peletier et al. (1995). For a 
nominal comoving galactic density of no = 0.002/if Mpc -3 
(eg. Efstathiou et al. 1988), these scale radii correspond to 
a range for r 9 (equation FlOI) : 0.01<r 9 <0.18. These ranges 
are consistent with those assumed in the intervening galaxy 
obscuration models of Heisler & Ostriker (1988) and Fall & 
Pei (1993). 

The values for (5, Zdust) were chosen to cover a range 
of evolution strengths for ro and tb respectively. To cover 
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a plausible range of dust formation epochs, we consider 
6 < Zdust < 20, consistent with a range of galaxy 'formation' 
epochs predicted by existing theories of structure formation 
(eg. Peebles 1989). The upper bound Zdust = 20 corresponds 
to the star formation epoch considered in the galaxy forma- 
tion models of Blain & Longair (1993b). 

We assume values for S similar to those implied by ob- 
servations of the space density of metal absorption systems 
from QSO spectra as a function of redshift (Sargent, Bok- 
senberg & Steidel 1988; Thomas & Webster 1990). These 
systems are thought to arise in gas associated with galaxies 
and their haloes and it is quite plausible that such systems 
also contain dust. Here we assume a direct proportionality 
between the amount of dust and heavy metal abundance in 
these systems. 

In general, evolution in the number of metal absorption 
line systems per unit z, that takes into account effects of 
cosmological expansion, can be parameterised: 



dN 
dz 



Ho 



n z Tvr (z) 2 (l + z)(l + 2q z)- 1/2 . 



(11) 



Evolution, such as a reduction in absorber numbers with red- 
shift, can be interpreted as either a decrease in the comoving 
number density n z , or effective cross-section Tvro(z) 2 . With 
our assumption of a constant comoving density n(z) = no, 
and an evolving dust scale radius ro as defined by equation 
(^), we have dN/dz oc (1 + z) 1 , where 7 = 0.5 + 2<5 for 
go = 0.5. Hence for no evolution, 7 = 0.5. 

Present estimates on the evolution of absorber numbers 
with redshift are poorly constrained. Thomas & Webster 
(1990) have combined several datasets increasing absorp- 
tion redshift ranges to give strong constraints on evolution 
models. For CIV absorption (AA1548, 155lA), which can be 
detected to redshifts z>3 in high resolution optical spec- 
tra, evolution has been confirmed for the highest equivalent 
width systems with WoJbO.GA. It is more likely that these 
systems are those associated with dust rather than the lower 
equivalent width (presumably less chemically enriched) sys- 
tems with Wo<0.3A which have a trend consistent with no 
evolution. Their value for the evolution parameter 7, for the 
highest equivalent width systems is — 0.1±0.5 at the 2a level. 
Converting this 2a range to our model parameter 5 using the 
discussion above, we assume the range: —0.5 < 8 < —0.05. 



3.3 Comparisons with QSO Absorption-Line 
Studies 

We can compare our assumed ranges in evolutionary pa- 
rameters: 6 < Zdust < 20 and —0.5 < S < —0.05 with recent 
determinations of the heavy element abundance in damped 
Ly-a absorption systems and the Ly-a forest to z ~ 3. The 
damped Ly-a systems are interpreted as the progenitors of 
galactic disks (Wolfe et al. 1986), and recent studies by Pet- 
tini et al. (1994; 1997) deduce metal abundances and dust- 
to-gas ratios at z ~ 1.8 — 2.2 that are ~ 10% of the local 
value. The Lyman forest systems however are more numer- 
ous, and usually correspond to gas columns > 10 7 times 
lower than those of damped Ly-a absorbers. High resolution 
metal-line observations by Songaila (1997) deduce metallic- 
ities <1.5% solar at z ~ 2.5 — 3.8. 

To relate these metallicity estimates to cosmic evolution 
in dust content as specified by our model, we must first note 
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Figure 2. Relative metallicity as a function of redshift. Regions 
within the solid and dashed curves represent the ranges predicted 
by our model for go = and 50 = 0.5 respectively. The filled and 
open data points with lcr error bars are mean observed values 
from Pettini et al. (1994) and Songaila (1997) respectively. 



that the metallicity at any redshift Z(z), is generally defined 
as the mass fraction of heavy metals relative to the total gas 
mass: Z(z) = Slm(z)/Q g (z). At all redshifts, we assume a 
constant dust-to-metals ratio, Qd{z)/Slm(z), where a fixed 
fraction of heavy elements is assumed to be condensed into 
dust grains. Therefore the metallicity Z(z), relative to the 
local solar value, Zq, can be written: 



z(z) _ Q d (z)n g (o) 



(12) 



z Q n d (o)n g (z)' 

From the formalism in section |2.l[ the mass density in dust 
relative to the local density, Q d {z) /Qd(0), can be determined 
and is found to be independent of the galaxy properties ro 
and tb, depending only on our evolution parameters, 5 and 
Zdust- This is given by 



n d {z) 



1 - 



ln(l + z) 



ln(l + Zdust) 



(1 + *) 



(13) 



The gas ratio, Q g (Q)/Q, g (z), is adopted from studies of the 
evolution in gas content of damped Ly-a systems. These 
systems are believed to account for at least 80% of the gas 
content in the form of neutral hydrogen at redshifts z>2 
(Lanzetta et al. 1991). We adopt the empirical fit of Lanzetta 
et al. (1995), who find that the observed evolution in Q g (z) 
is well represented by fl g (z) = S7 g (0) exp(az), where a = 
0.6L0.15 and 0.83L0.15 for qo = and go = 0.5 respectively. 

Figure || shows the range in relative metallicity implied 
by our evolutionary dust model (equations [[2] and |l|) as 
a function of redshift for two values of go- The solid and 
dashed lines correspond to respectively go = and go = 0.5 
and the regions between these lines correspond to the ranges 
assumed for our assumed model parameters: 6 < Zdust < 20 
and —0.5 < S < —0.05. For comparison, the mean metal- 
licities Z ~ 0.1 Zq and Z ~ 0. 01 Zq observed in damped 
Ly-a systems at z ~ 2.2 and the Lyman forest at z>2.5 re- 
spectively are also shown. These agree well with our model 
predictions, suggesting that our model assumptions will pro- 
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Figure 3. Optical depth probability distribution functions p(r\z) to redshifts z =1, 3 and 5, where r is the total optical d ept h observed 
in the S-band. Two different sets of galaxy parameters (r g , tb) are considered: (a) (0.2,4) and (b) (0.01,0.5) (see section p,l| ). For each 
of these, we show four evolutionary models specified by (5, z dust ). 'No evolution' corresponds to S = and z dust = oo and the 'Strongest 
evolution' to 8 = —0.5 and z dust = 6. 



vide a reliable measure of dust evolution which are at least 
compatible with other indirect estimates. 



4 RESULTS AND ANALYSIS 

Using the formalism of Masci & Webster (1995) and replac- 
ing the parameters tb and ro by their assumed redshift de- 
pendence as denned in section 2.1, Fig. ^| shows probability 
density functions p(r | z) for the total optical depth up to 
redshifts z =1, 3 and 5. Results are shown for two sets of 
galaxy parameters (r g , tb), with four sets of evolutionary 
parameters (8, Zdust) for each. 

The area under any normalised curve in Fig. |§] gives the 
fraction of lines-of-sight to that redshift which have optical 
depths within some interval — > r max . Towards high red- 
shifts, we find that obscuration depends most sensitively on 
the parameter r B , in other words, on the covering factor of 
absorbers (equation . Figure ^| shows that as the amount 
of dust at high redshift decreases, ie., as S and Zdust de- 
crease, the curves show little horizontal shift towards larger 
optical depths from z — 1 to z = 5. A significant shift be- 
comes noticeable however for the weaker evolution cases, 



and is largest for 'no evolution' (solid lines). This behaviour 
is further investigated below. 

In order to give a clearer comparison between the 
amount of obscuration and strength of evolution implied by 
our model parameters (r 9 , tb, 8, z dust ), we have calculated 
the mean and variance in total optical depth as a function 
of redshift. Formal derivations of these quantities are given 
in the appendix. Here we briefly discuss their general depen- 
dence on the model parameters. 

A quantity first worth considering is the number of 
galaxies intercepted along the line-of-sight. In a go = 0.5 
(A = 0) universe, the average number of intersections within 
a scale length ro of a galaxy's center by a light ray to some 
redshift is given by 



1.5 + 25 



-1 



(14) 



Where S and r g are defined in equations (|^) and ( pTj| ) respec- 
tively. 

In the case where we have no-evolutton, ie. where 8 — 
and z ma x = oo, and for a dust law that scales inversely with 
wavelength (ie. £a oc 1/A which is a good approximation at 
A>2500A), exact expressions follow for the mean and vari- 
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Figure 4. Behaviour in mean reddening, (t), as a function of redshift for a range of model parameters (r„, Tg) and (5, z dust ). (a) For 
(r 9 , Tg)=(0.2,4) and (0.2,0.5), (b) Same as (a) but for t 9 = 0.01, (c) Redshift dependence of mean reddening in no-evolution model for 
a range of parameters (t 9 , Tg). (d) Scaling of the mean reddening with respect to the evolutionary parameters (5, z dust ) with (t 9 , Tg) 
fixed at (0.2,4). 



ance in total optical depth along the line-of-sight. The mean 
optical depth can be written: 

f(z) = 0.8 r g tb [(1 + z) 2 ' 5 - l] , (15) 
and the variance: 

a*(z) = 0.57 t 9 tb 2 [(1 + z) 3 - 5 - l] . (16) 

The variance (equation ^) or 'scatter' about the mean 
to some redshift provides a more convenient measure of red- 
dening. The mean optical depth has a simple linear depen- 
dence on the parameters t 9 and tb and thus gives no in- 
dication of the degree to which each of these parameters 
contributes to the scatter. As seen from the probability dis- 
tributions in Fig. ^ there is a relatively large scatter about 
the mean optical depth to any redshift. From equation (|l6|), 
it is seen that the strongest dependence of the variance is on 
the central absorber optical depth tb- Thus, larger values of 
tb (which imply 'harder-edged' disks), are expected to intro- 
duce considerable scatter amongst random individual lines 
of sight, even to relatively low redshift. 

In Fig. |^, we show how the mean optical depth varies 
as a function of redshift for a range of evolutionary pa- 
rameters. 'Strong evolution' is characterised by S = —0.5, 



Zdust = 6 (dot-dashed curves), as compared to the 'no', 
'weak' and 'moderate' evolution cases indicated. The mean 
optical depth flattens out considerably towards high redshift 
in the strong evolution case, and gradually steepens as 5 and 
Zdust are increased. Note that no such flattening is expected 
in mean reddening for the no evolution case (Fig. gj). The 
mean optical depth to redshifts «<;1 in evolution models can 
be reduced by factors of at least three, even for low to mod- 
erately low evolution strengths. 

Figure ^d shows the scaling of mean optical depth with 
respect to the evolutionary parameters. It is seen that red- 
dening depends most sensitively on the parameter 8, which 
controls the rate of evolution in galactic dust scale radius 
ro. A similar trend is followed in Fig. ^, which shows the 
dependence of variance in optical depth on evolution as a 
function of redshift, for fixed (r 9 , tb)- Considerable scat- 
ter is expected if the dust radius of a typical galaxy evolves 
slowly with cosmic time as shown for the 'weakest' evolution 
case 8 = —0.05 in Fig. ^. 

Our main conclusion is that the inclusion of evolution 
in dust content, by amounts consistent with other indirect 
studies can dramatically reduce the redshift dependence of 
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Figure 5. Variance (a^) in optical depth as a function of red- 
shift showing scaling with respect to the evolution parameters (<5, 
z dust)- ( T gi T s) are fixed at (0.2,4). (a) <5 fixed at -0.05 and Zdust 
is varied, (b) z^ ugt fixed at 10 and S is varied. 



total reddening along the line-of-sight to z>l, contrary to 
non-evolving models. 



5 IMPLICATIONS ON QSO NUMBER 
COUNTS 

There are numerous observations suggesting that the space 
density of bright quasars declines beyond z ~ 3 (Sandage 
1972; Schmidt, Schneider & Gunn 1988). This has been 
strongly confirmed from various luminosity function (LF) 
estimates to 2 ~ 4.5 (Hartwick & Schade 1990; Pei 1995 
and references therein), where the space density is seen to 
decline by at least an order of magnitude from z — 3 to 
z — 4. Heisler & Ostriker (1988) speculate that the decline 
may be due to obscuration by intervening dust, which re- 
duces the number of quasars observed by ever-increasing 
amounts towards high z. The results of Fall & Pei (1993) 
however show that the observed turnover at z ~ 2.5 and de- 
cline thereafter may still exist once the effects of intervening 
dust (mainly associated with damped Lya systems) are cor- 
rected for. Since no evolution in dust content was assumed 
in either of these studies, we shall further explore the effects 
of intervening dust on inferred quasar evolution using our 
evolutionary galactic dust model. 

Since we are mainly interested in "bright" quasars 



< 



26) at high redshifts, a single power-law for the 



observed LF should suffice: 
(f>o{L,z) = (j>, (z)L~ p ~ 1 , 



(17) 



where f3 ~ 2.5. This power law model immensely simplifies 
the relation between observed and "true" LFs (corrected for 
obscuration by dust). In the presence of dust obscuration, 



inferred luminosities will be decreased by a factor of e _T . 
Since there is a probability p(r | z) of encountering an optical 
depth r as specified by our model (see Fig. ^|), the observed 
LF can be written in terms of the true LF, d>± as follows: 



4>o(L, z) 



dr (f> t (e T L, z)e T p(r \ z) 



(18) 



The extra factor of e T in equation (jl8|) accounts for a de- 
crease in luminosity interval dL in the presence of dust. 
Equations (0) and (jyj) imply that the true LF can be writ- 
ten 

<t>t{L,z) = <f>, t (z)L- -\ (19) 
and the ratio of observed to true LF normalisation as 

4>*o(z) _ f°° j_„-/3t 



dr e p(r I z). 



(20) 



The observed comoving density of quasars brighter than 
some absolute magnitude limit Mu m as a function of redshift 
is computed by integrating the LF: 



N (z\Mb < M Hm ) = / dL 4>o(L, z). (21) 

J L lim, = L(M lim ) 

Thus, the true comoving number density jV t ,can be eas- 
ily calculated by replacing (f> in equation (Elh by cf>t = 
(4>*t/4>*o)4'o leading to the simple result: 



N t (z\M B < M lim ) 



<P*o(z) 
<t>*M 



N (z\M B <M Um ), (22) 



where the normalisation ratio is defined by equation (|20|), 

Figure ^| shows both the observed and true comoving 
density of bright quasars (with Mb < —26) as a function 
of redshift. The observed trends are empirical fits deduced 
by Pei (1995). The true comoving density in all cases was 
determined by assuming relatively 'weak' evolution in the 
dust properties of intervening galaxies. Two sets of galactic 
dust parameters for each qo defined by (rs,ro) = (1, lOkpc) 
(Figs a and c) and (tb,to) = (3, 30kpc) (Figs b and d) 
are assumed. We shall refer to these as our "minimal" and 
"maximal" dust model respectively which bracket the range 
of parameters observed for local galaxies. 

Comparing the 'true' QSO redshift distribution with 
that observed, two features are apparent. First, the true 
number density vs. z relation has qualitatively the same be- 
haviour as that observed. No flattening or increase in true 
quasar numbers with z is apparent. Second, there appears 
to be a shift in the redshift, z pea k, where the quasar density 
peaks. This shift is greatest for our maximal dust model 
where z pea k is increased by a factor of almost 1.5 relative 
to that observed. This implies that the bulk of quasars may 
have formed at earlier epochs than previously inferred from 
direct observation. 

Our predictions for QSO evolution, corrected for ob- 
scuration by 'evolving' intervening dust differs enormously 
from that predicted by Heisler & Ostriker (1988). The major 
difference is that these authors neglected evolution in dust 
content with z. As shown in Fig. ^, non-evolving models lead 
to a rapid increase in dust optical depth with z and hence 
this will explain their claim of a continuous increase in the 
true QSO space density at z > 3. As shown in Fig. ^, the 
inclusion of even a low-to-moderately low amount of evolu- 
tion in dust content dramatically reduces the excess number 
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Figure 6. Comoving number density of quasars with Mg < —26 
as a function of redshift. Observed trends (solid curves) are taken 
from the empirical fits of Pei (1995) while dashed curves corrects 
these trends for obscuration by dust. These are predicted assum- 
ing our evolving intervening galactic dust model with Tg = 1 and 
ro = lOkpc (Figs a and c) and tb = 3 and ro = 30kpc (Figs b 
and d). In all cases, we have assumed relatively "weak" evolution 
in dust content with z, defined by the parameters: z^ ust = 20 and 
<5 = -0.05. 



of quasars at z > 3 than predicted by Heisler & Ostriker 
(1988). 

We find that there is no significant difference in the 
characteristic timescale, tQso for QSO formation at z > 
Zpeak, where 



, N 

tQSO — — 



1.5 Gyr, 



(23) 



is found for both the observed and dust corrected results 
in Fig. |(| We conclude that the decline in space density of 
bright QSOs at redshifts z > 3.5 is most likely to be real 
and an artifact of an intrinsic rapid turn-on of the QSO 
population with time. This is consistent with estimates of 
evolution inferred from radio-quasar surveys where no bias 
from dust obscuration is expected (eg. Dunlop & Peacock 
1990). 

An increased space density of quasars at redshifts z > 3 
predicted by correcting for dust obscuration has implica- 
tions for theories of structure formation in the Universe. Our 
minimal dust model (Figs. |^a and c) predicts that the true 
space density can be greater by almost two orders of mag- 
nitude than that observed, while our maximal dust model 
(Figs. [j]b and d) predicts this factor to be greater than 5 



a 


- b ^----^ 

/ \- 


/ "observed" \\ - 

"true" \\ 




q = 
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orders of magnitude. These predictions can be reconciled 
with the quasar number densities predicted from hierarchi- 
cal galaxy formation simulations involving cold-dark matter 
(eg. Katz et al. 1994). It is found that there are > 10 3 times 
potential quasar sites at z > 4 (associated with high den- 
sity peaks) than required from current observations. Such 
numbers can be easily accommodated by our predictions if 
a significant quantity of line-of-sight dust is present. 

To summarise, we have shown that with the inclusion of 
even weak to moderately weak amounts of evolution in dust 
content with z, the bias due to dust obscuration will not 
be enough to flatten the true redshift distribution of bright 
quasars beyond z = 3. A significant excess however (over 
that observed) in quasar numbers is still predicted. 



6 DISCUSSION 

Our model predictions may critically depend on the dust 
properties of individual galaxies and their assumed evolu- 
tion. For instance, is it reasonable to give galaxies an ex- 
ponential dust distribution? Such a distribution is expected 
to give a dust covering factor to some redshift considerably 
larger than if a clumpy distribution were assumed (Wright, 
1986). A clumpy dust distribution (for spirals in particular) 
is expected, as dust is known to primarily form in dense, 
molecular star- forming clouds (Wang 1991 and references 
therein). 

As noted by Wright (1986), "cloudy disks" with dust in 
optically-thick clumps can reduce the effective cross-section 
for dust absorption by at least a factor of five and hence, are 
less efficient at both obscuring and reddening background 
sources at high redshift. A dependence of the degree of dust 
'dumpiness' on redshift, such as dust which is more dif- 
fuse at early epochs and becomes more clumpy with cosmic 
time is unlikely to affect the results of this paper. This will 
only reduce the effective cross-section for absorption to low 
redshifts, leaving the effects to high redshift essentially un- 
changed. The numbers of reddened and/or obscured sources 
at high redshift relative to those expected in non-evolving 
dust models however will always be reduced, regardless of 
the dependence of aborption cross-section on redshift. 

Observations of the optical reddening distribution of 
quasars as a function of redshift may be used to test our 
predictions. Large and complete radio-selected samples with 
a high identification rate extending to high redshifts how- 
ever are required. The reason for this is that first, radio 
wavelegths are guaranteed to have no bias against obscura- 
tion by dust, and second, the statistics at high redshift need 
to be reasonably high in order to provide sufficient sampling 
of an unbiased number of random sight-lines. 

The sample of Drinkwater et al. (1997) contains the 
highest quasar fraction (Jj70%) than any existing radio sam- 
ple with a redshift distribution extending to z ~ 4. A large 
fraction of sources appear very red in B — K colour com- 
pared to quasars selected by optical means. The dependence 
of B — K colour on redshift is relatively flat which may at 
first appear consistent with the predictions of figure ^| al- 
though the fraction of sources identified with z>2 is only 
~ 5%. Also, this sample is known to contain large numbers 
of sources which are reddened by mechanisms other than 
dust in the line-of-sight (eg. Serjent & Rawlings 1996). The 
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role of dust, reddening the optical-to-near-IR continua of 
radio-selected quasars, and whether it is extrinsic or not still 
remains a controversial issue. One needs to isolate the intrin- 
sic source properties before attributing any excess reddening 
to line-of-sight dust. Optical follow-up of sensitive radio sur- 
veys that detect large numbers of high redshift sources with 
known intrinsic spectral properties will be necessary to reli- 
ably constrain the rate of evolution in cosmic dust. 



7 SUMMARY AND CONCLUSIONS 

In this paper, we have modelled the optical depth in galactic 
dust along the line-of-sight as a function of redshift assum- 
ing evolution in dust content. Our model depends on four 
parameters which specify the dust properties of local galax- 
ies and their evolution: the exponential dust scale radius 
ro, central B-band optical depth tb, "evolution strength" 5 
where ro(z) — ro(l + z) s , and Zdust - a hypothesised dust 
formation epoch. Our evolution model is based on previous 
studies of the formation of heavy metals in the cold dark 
matter scenario of galaxy formation. 
Our main results are: 

1. For evolutionary parameters consistent with exist- 
ing studies of the evolution of metallicity deduced from 
QSO absorption-line systems, a significant "flattening" in 
the mean and variance of observed B-band optical depth 
to redshifts z > 1 is expected. The mean optical depth to 
z~>l is smaller by at least a factor of 3 compared to non- 
evolving model predictions. Obscuration by dust is not as 
severe as shown in previous studies if effects of evolution are 
accounted for. 

2. By allowing for even moderately low amounts of evo- 
lution, line-of-sight dust is not expected to significantly af- 
fect existing optical studies of QSO evolution. Correcting for 
dust obscuration, evolving dust models predict the 'true' 
(intrinsic) space density of bright quasars to decrease be- 
yond z ~ 2.5 as observed, contrary to previous non-evolving 
dust models where a continuous monotonic increase was pre- 
dicted. 

3. For moderate amounts of evolution, our models pre- 
dict a mean observed B-band optical depth that scales as 
a function of redshift as f a (1 + z) - 1 . For comparison, 
evolving models predict a dependence: f oc (1 + z) 2 ' 5 . We 
believe future radio surveys of high sensitivity that reveal 
large numbers of optically reddened sources at high redshift 
will provide the necessary data to constrain these models. 
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APPENDIX A: DERIVATION OF MEAN 
OPTICAL DEPTH 

Here we derive expressions for the mean and variance in total 
optical depth as a function of redshift in ou r evolutionary 
galactic dust model discussed in section 3.2. The galaxies 



are modelled as exponential dusty disks, randomly inclined 
to the line-of-sight. 

We first derive the average number of galaxies inter- 
cepted by a light ray emitted from some redshift z (ie. equa- 
tion |l^). Given a 'proper' number density of galaxies at some 
redshift n g (z), with each galaxy having an effective cross- 
sectional area \xo as viewed by an observer (u is a random 
inclination factor, where /i = cos 9 and 8 is the angle be- 
tween the sky plane and the plane of a galactic disk), the 
average number of intersections of a light ray along some 
path length ds will be given by 



dN = n g (z) fj,a ds. 



(Al) 



In an expanding universe we have n g = noil + z) 3 , where 
no is a local comoving number density and is assumed to be 
constant. Units of proper length and redshift are related by 



ds 
dz 



\H ) (1 



(l + z) 2 (l + 2 (?0 z) 1 /2 



(A2) 



(Weinberg 1972). The effective cross-section projected to- 
wards an observer for a randomly inclined disk is found by 
averaging over the random inclination factor fi, where \i is 
randomly distributed between and 1, and integrating over 
the exponential profile assumed for each disk with scale ra- 
dius ro(z) (see equations [j] and|^). The product [io in equa- 
tion (Al) is thus replaced by 



/i dfi 



(A3) 



Thus from equation (Al), the average number of intersec- 



tions to some redshift z is given by 

ds 



N{z) = 



fia 



^ {-£) dz ' 

(l + z') 1+2S 



c 



(l + 2g z') 1/i 



dz' . 



(A4) 



With T q defined by riQ-Kr\ (jf-J, this directly leads to equa- 
tion (Q) for q = 0.5. 

The mean optical depth f is derived by a similar argu- 
ment. If to (z) is the optical depth observed through a face 
on galaxy at some redshift z (equation ^), then a galactic 
disk inclined by some factor [s, will have its optical depth 
increased to To{z)/fi. Multiplying this quantity by equation 
(Al), the extinction suffered by a light ray along a path 
length ds is given by 



dr = rig(z) a tq(z) ds. 



(A5) 



Thus the mean optical depth to some redshift z can be cal- 
culated from 



f(z) = 



an g (z')T {z) {j^j dz' 



(A6) 



Giv en n g (z), (^f ) and a (from the integral over r in equation 
A3) above, and to(jz') from equation (fef), the mean optical 



depth follows the general form 

(1 + 2') 1 + 2 " 



f(z) 



2r » TB J (i + 2 qo z>y/^{—) 



ln(l + z') 



ln(l + Zdus 



dz' . 



(A7) 



Similarly, the variance in the optical depth distribution is 
defined as follows: 



4(Z) = (r 2 > - (r) 2 = 



,(*W) dz'. (A8) 



In terms of our model dependent parameters, this becomes 
' (l + z') 1+2S A £ 



o 2 {z) = 2t 9 t b 2 



(1 + z)— 2 ( Xb \ 
(1 + 2 (?0 2') 1/2 VI + zj 



ln(l + z') 



ln(l + Zdust) 



dz' . 



(A9) 



